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The use of chiral perturbation theory in extracting physical K — > %n matrix elements from matrix 
elements calculated at unphysical kinematics is outlined. In particular, the possibility of utilizing 
pions with non-zero momentum in the final state, and of using partial quenching is discussed. 
Preliminary (not physically normalized) AI = 3/2 (27,1) K — > nn matrix elements are calculated 
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fermions, with an inverse lattice spacing of a -1 = 1.729(28) GeV. Effective mass plots are pre- 
sented for a light sea quark mass of m^ ea = 0.005, and various valence quark masses. The plateaux 
are fit and E Kn — trig is extracted. 
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1. Introduction 



The decay K ^ KK is important to simulate on the lattice because QCD effects come into play 
since the typical energies are smaller than the scale Aqcd- In particular, the direct CP violating 
parameter Re(e'/e) can be found from calculations of K — * %% matrix elements (dJ, ^J. In order 
to obtain reasonable precision in these calculations we plan to use 2+1 flavors of domain wall 
fermions (DWF) on a 24 3 x 64, L, = 16 lattice. 

Accurate simulation of physical pions would require a larger box size than is currently avail- 
able, therefore we make use of the tools of chiral perturbation theory (^PT). ^PT allows one to 
make predictions of the dependence of matrix elements on quark masses when they are close to the 
chiral limit. It can even be applied to the case of unphysical kinematics. Thus the best strategy is 
to simulate K — > %tz at unphysical kinematics, fit the results to jYY, and then use #PT to make a 
prediction for the physical values of the matrix elements. 



2. Four Quark Operators and #PT 

The weak interactions are included in the lattice QCD simulation by evaluating matrix ele- 
ments of the effective Hamiltonian [|3|, ^] 

^S=X = ^Z V CKMCi{pi)Qi (2-1) 

where c,-(ju) are Wilson coefficients and {Qi,i = 1, 10} are four quark operators. Therefore we 
are interested in calculating matrix elements of the four quark operators <2, between a K and a 

3/2 1 /2 

Tin state. These operators can be split into AI = 3/2 and AI = 1/2 parts, called Q t and Q/ 
respectively, where AI is the change in isospin induced by the operator. They can then be further 
classified by how they transform under the chiral SU (3)l x SU (3)r symmetry, and the transforma- 
tion properties (27,1), (8,8), and (8,1) are all found among various of the operators ^ 
In %PT an effective Lagrangian [§] is written in terms of the field 



£ = exp 



2i(j) a X a 
f 



(2.2) 



where the 0" are the real pseudo-scalar meson fields. The leading order part of the effective La- 
grangian is 

%w = ^Tr[d^n + i^Tr[x^ + ^ X ] (2.3) 

™2 m 2 



where % = diag{m u ,m d ,m s ) and B = — *±- = = To represent the four quark oper- 



nr 
K 

m u +m d m„+m s m d +m s 

ators in %PT, one looks at the operators that have definite transformation properties under isospin 
and SU (3)l x SU (3)#, and forms operators out of the £ field that transform in the same way. In 
general it will be possible to form more than one such operator, and thus a linear combination of 
all of these operators in which each operator is multiplied by an arbitrary coefficient is taken to 
represent the four quark operator [0]. These arbitrary coefficients are called low energy constants 
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(LECs). A %PT expression for a K — > %% matrix element of a four quark operator at next to lead- 
ing order (NLO) and with pions in the final state having zero three-momentum, generally contains 
polynomials in the meson masses with order m 2 and order m 4 terms, as well as non-analytic terms 
in the meson masses known as chiral logarithms. The coefficients of these terms contain LECs and 
other parameters in the Lagrangian. 

3. Extraction of the LECs 

Matrix elements of four quark operators can be computed on the lattice for several different 
quark masses. The results can be fit to the %PT formulae and the fit will yield LECs. The strategy 
is to compute these matrix elements at unphysical kinematics in order to extract the LECs that 
appear in the %PT expression for the physical matrix elements. However, in order to be able to 
determine the necessary LECs uniquely with a limited number of gauge field ensembles one must 
either consider pions with non-zero momenta, or one must resort to partial quenching in which the 
masses of the quarks in the fermion determinant (sea quarks) are different from the masses of the 
propagating quarks (valence quarks). 

3.1 Pions with Non-Zero Momenta 

%PT formulae for K —> %% matrix elements with pions having non-zero momentum have been 
worked out by Lin et. al. [0] and Laiho and Soni [J|J|. One drawback of this method is that data with 
non-zero momentum tends to be very noisy, especially as the momentum increases. There exist 
methods for dealing with this such as antiperiodic, and in general twisted boundary conditions. 

3.2 Partial Quenching 

The largest 24 3 x 64, L 4 = 16 RBC/UKQCD 2+1 flavour dynamical lattice ensembles currently 
available have sea quark masses m^ ea = m^ ea = {0.005,0.01} and mf a = 0.04. In general it is 
necessary to vary both the light quark and strange quark valence masses in order to extract LECs. 
Laiho and Soni [||] have treated the case of partially quenched %PT at NLO with sea quarks of equal 
mass. However, partially quenched %PT formulae at NLO for the case of unequal sea quark masses 
do not yet exist in the literature and are in the process of being calculated by Christopher Aubin, 
Shu Li, and Jack Laiho. The current plan is to focus primarily on the partial quenching technique. 
Non-zero momenta can then be incorporated as a consistency check of the LECs obtained from the 
former method. 

Table 1 gives a list of sets of quark masses that would be sufficient to extract the necessary 
LECs from an analysis of [gj. However, it is not guaranteed that these masses, with the sea strange 
quark mass changed to mf 3 ' = 0.04 in all cases, will still be sufficient when one considers %PT 
with unequal sea quark masses. This analysis was done merely as an exercise to get a feeling for 
the number of combinations of masses that would be needed. 
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Table 1: A list of sets of quark masses at which to evaluate several matrix elements, from which the LECs 
needed to calculate the physical K — > Tin matrix elements can be extracted. 



mf a = mf a 


m™ 1 


< 


0.01 


0.01 


0.01 


0.02 


0.02 


0.02 


0.02 


0.01 


0.01 


0.005 


0.035 


0.005 


0.01 


0.02 


0.01 


0.02 


0.03 


0.01 



4. Matrix Element Calculations 

So far K — > %% matrix elements of the AI = 3/2 (27,1) four quark operator have been calcu- 
lated. The operators that have a non-trivial AI = 3/2 part that transforms as a (27,1) are Q\, Q2, 

3/2 3/2 3/2 3/2 

Q9, Qw- It turns out by Fierz symmetry that £2/ , Q 2 , Qg » 610 are all proportional to a single 
operator Q&MW) nTjJ 

0(27,0(3/2) = 3Q 3/2 = 2Q 3/2 = ^(1 _ f) da u b f(\ - ^)u b +f«y M (l - tV^I " T^K 

-^7m(i - ^KAr^i - (4-1) 

Furthermore, the Wigner-Eckhart theorem can be used to write 

{% + % + \^ 21 ^l 2 \K + ) = -^(n + n°\Q^- l ^ 2 \K + ) (4.2) 

where 

0/(27,1)^/2 = fflyM(1 _ yS^g^j _ ^ (4 3) 

The matrix element on the left hand side of Equation [O] is easier to deal with, and from this matrix 
element one can see that only one diagram contributes to AI = 3/2 (27,1) matrix elements. This 
diagram is shown in Figure [j] 



K+ 
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Figure 1: The one diagram that contributes to the AI = 3/2 (27,1) K — > Tin matrix element. 

The AI = 3/2 (27,1) matrix elements have been calculated on the RBC/UKQCD 24 3 x 64, 
L s = 16 lattices, using 2+1 dynamical flavors and domain wall fermions, with an inverse lattice 
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Table 2: Values of E Kn — m^ obtained by fitting the plateaux of the effective mass plots. Errors obtained 
from the jackknife method are quoted. Here mf* = 0.005 and mf a — 0.04. 



m s 


mi 


- m K 


m s 


mi 




m s 


mi 


E n % - m K 


0.04 


0.04 


0.4277(9) 


0.03 


0.01 


0.1573(20) 


0.005 


0.005 


0.1835(25) 


0.04 


0.03 


0.3459(9) 


0.02 


0.01 


0.1930(20) 


0.04 


0.001 


-0.0483(46) 


0.03 


0.03 


0.3737(10) 


0.01 


0.01 


0.2333(20) 


0.03 


0.001 


-0.0135(39) 


0.04 


0.02 


0.2496(13) 


0.04 


0.005 


0.0420(28) 


0.02 


0.001 


0.0275(38) 


0.03 


0.02 


0.2795(14) 


0.03 


0.005 


0.0758(26) 


0.01 


0.001 


0.0768(35) 


0.02 


0.02 


0.3115(14) 


0.02 


0.005 


0.1138(25) 


0.005 


0.001 


0.1076(30) 


0.04 


0.01 


0.1249(21) 


0.01 


0.005 


0.1579(25) 


0.001 


0.001 


0.1371(29) 



spacing of a- 1 = 1.729(28) GeV. (These lattices are described in more detail in [|TT|]). Wall sources 
at t = 5 and t = 59 are used for the kaon and two pions respectively. Averages over periodic and 
antiperiodic boundary conditions are performed in order to double the effective time length and 
prevent contamination by 'around the world' propagation. Matrix elements are evaluated with 
the four quark operator at different times between t = 5 and ? = 59. Calculations are done with 
sea quark masses mf a = 0.04, mf a = m s J a = 0.005,0.01, and valence quark masses in the set 
{0.001,0.005,0.01,0.02,0.03,0.04} and all possible combinations such that mf > mf x (where 
mi = m u = m d ). 

Figures |2| - |] show effective mass plots of the matrix element as a function of the time at 
which the four quark operator is located. In these plots the sea quark masses are mf a = 0.04 and 
mf a = 0.005, and the valence light quark mass is held fixed while the valence strange quark mass 
is varied. The plateaux are fitted in a range t m [ n <t< f max which is different for each curve, and 
the fits are indicated by bold lines in the plots. 

The effective mass is defined as 

i^Ur) (44) 

where C(?) is the value of the matrix element at time t. For times far away from the sources and 
sinks we expect that 

C{t)^Aexp[(E K7l -m K )t] (4.5) 

and thus that 

m ef f(f) « E KK - m K (4.6) 

Therefore we expect the value of the plateau to be E %% — m^. Results for this quantity obtained by 
fitting the plateaux are given in Table 2. 

5. Conclusion and Future Plans 

Calculations of AI = 3/2 (27 , 1 ) K — > %% matrix elements have been done on the RBC/UKQCD 
2+1 flavour dynamical 24 3 x 64, L s = 16 lattices. To obtain physically normalized matrix elements 
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Figure 2: Effective mass plots of the K — > 7T7Z; matrix element as a function of the time at which the four 
quark operator is located. The sea quark masses are mf a = 0.005 and mf a — 0.04 in all plots here. The 
valence quark masses are m™ 1 = 0.04 (left), w™ 1 = 0.03 (right), and various values of m™ 1 are shown as 
different curves. A horizontal line (bold) is fitted in a range t m in — t — ^max which is different for each curve. 




Figure 3: Effective mass plots for valence quark masses mj a = 0.02 (left), m™ 1 = 0.01 (right), and various 
values of m™ 1 shown as different curves. 
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it is necessary to use kaon, pion, and two pion correlators, the latter of which are in the process of 

3/2 3/2 

being calculated. Then it will be possible to move onto matrix elements of Q 7 and <2 8 which 
transform as (8,8). The next step is to calculate matrix elements of the Al = 1/2 operators, however 
these present additional computational difficulties related to large vacuum subtractions. 

These matrix elements on the lattice must also be renormalized to their continuum values. 
This can be done using the technique of non-perturbative renormalization (NPR). Work on NPR 
for K — ► KK has been done by Shu Li. 

Finally, the matrix elements must be fit to NLO %PT expressions as soon as they become 
available. This should yield the LECs necessary to calculate the physical value of the K — > KK 
matrix element. We would also like to calculate some matrix elements with pions that have non- 
zero momenta, for example the first non-vanishing value of momentum allowed by the lattice, and 
use this as a consistency check of the results obtained from partial quenching. The RBC/UKQCD 
is also working on another large lattice of dimensions 32 3 x 64. Other lattice sizes could be used 
to perform a continuum extrapolation. 
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